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Abstract 

While the role of estrogen receptor-related receptor alpha (ERRa) in chondrogenesis has been investigated, the involvement 
of ERR gamma (ERRy) has not been determined. To assess the effect of increased ERRy activity on cartilage development 
in vivo, we generated two transgenic (Tg) lines overexpressing ERRy2 via a chondrocyte-specific promoter; the two lines 
exhibited ~3 and ~5 fold increased ERRy2 protein expression respectively in El 4.5 Tg versus wild type (WT) limbs. On 
postnatal day seven (P7), we observed a 4-10% reduction in the size of the craniofacial, axial and appendicular skeletons in 
Tg versus WT mice. The reduction in bone length was already present at birth and did not appear to involve bones that are 
derived via intramembranous bone formation as the bones of the calvaria, clavicle, and the mandible developed normally. 
Histological analysis of P7 growth plates revealed a reduction in the length of the Tg versus WT growth plate, the majority 
of which was attributable to a reduced proliferative zone. The reduced proliferative zone paralleled a decrease in the 
number of Ki67-positive proliferating cells, with no significant change in apoptosis, and was accompanied by large cell-free 
swaths of cartilage matrix, which extended through multiple zones of the growth plate. Using a bioinformatics approach, 
we identified known chondrogenesis-associated genes with at least one predicted ERR binding site in their proximal 
promoters, as well as cell cycle regulators known to be regulated by ERRy. Of the genes identified, Col2al, Agg, Pthlr, and 
Cdknlb (p27) were significantly upregulated, suggesting that ERRy2 negatively regulates chondrocyte proliferation and 
positively regulates matrix synthesis to coordinate growth plate height and organization. 
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Introduction 

The bones of the axial and appendicular skeleton arise from 
condensations of chondrogenic cells that lay down a cartilaginous 
scaffold, which is then remodeled to give rise to the ossified bone. 
The longitudinal growth of endochondral bones is driven by 
continued chondrogenesis in the growth plate, which can be 
divided along the longitudinal axis of the bone into distinct zones 
comprising resting, proliferating, and post-mitotic chondrocytes, 
respectively, from the articular surface. The process of chondro- 
genesis is a highly orchestrated proliferation-differentiation 
sequence that is regulated by a number of signaling pathways 
and feedback loops. For example, the homeobox transcription 
factor, SRY-related high-mobility-group box 9 (SOX9), is the 
primary determinant of chondrogenesis and is required for the 
initial commitment of mesenchymal stem cells to the chondrogenic 
lineage [1]. The hedgehog protein family member, Indian 
hedgehog (IHH), is also crucial to chondrocyte proliferation and 
differentiation. It is expressed by prehypertrophic cells and binds 



to its receptor Patched- 1 (PTC-1), which, in turn, activates 
signaling pathways to promote chondrocyte proliferation [2] . IHH 
also forms a negative regulatory feedback loop with parathyroid 
hormone-related protein (PTHRP) to delay chondrocyte hyper- 
trophy and increase the pool of proliferating chondrocytes [3]. On 
the other hand, the transcription factor RUNX2 promotes the 
differentiation of chondrocytes from proliferation to hypertrophy 
[4]. 

Certain transcription factors belonging to the nuclear hormone 
receptor family are also involved in chondrocyte differentiation. 
These include the two estrogen-binding receptors, estrogen 
receptor alpha and beta (ERot (NR3A1) and ERf3 (NR3A2) 
respectively), and recent reports highlight a role for ERot in the 
fusion or slowing down of growth plate chondrogenesis at puberty 
in humans and mice. For example, in a cartilage-specific ERot- 
deleted mouse, appendicular bones developed normally, but 
exposure to high levels of estrogen failed to reduce bone length 
as it did in wild type (WT) mice, indicating that ERot was required 
for the natural deceleration of bone growth that occurs in mice 
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Figure 1. Col2::ERRy2 transgenic mouse generation and protein expression analysis. (A) ERR and ER gene expression (relative to L32) in 
adult mouse cartilage. Graphs represent mean ± SD from a minimum of 3 independent samples. (B) Schematic of the transgene used to express the 
long protein isoform, ERRy2, in cartilage. Whole mount lacZ staining of a WT and transgenic 14.5 dpc embryo with a close-up of the forelimbs 
showing robust staining in the developing cartilage. Some nonspecific staining along the neural tube is evident. (C) Western blots and resulting 
quantification showing the levels of ERRy and p-ACTIN proteins in 14.5 dpc limbs of the two transgenic lines compared to WT. Graphs represent the 
mean ± SD from at least three independent Western blots. **p<0.01. 
doi:1 0.1 371 /journal.pone.0081 51 1 .g001 



upon sexual maturity [5]. Conversely, a mouse line that expressed 
a constitutively active form of ERa in cartilage exhibited fewer 
proliferating cells in the growth plate and reduced bone length [6] . 

Three orphan nuclear receptor genes related to the ERs 
comprise the estrogen receptor-related receptor (ERR) family: 
alpha, beta and gamma (NR3B1, NR3B2, and NR3B3, respec- 
tively) [7] . These genes share a high degree of similarity with the 
ERs, including 67% identity in the DNA-binding domain (DBD), 
but are unable to bind estrogen [8], With their similarity in their 
DBD, it is not surprising that there is considerable cross-talk at the 
level of gene regulation between the ERs and the ERRs. However, 
X-ray crystallography studies have clearly shown that, unlike the 
ERs, the ERRs assume an active state without a ligand bound to 
the ligand binding domain (LBD) [9,10]. Consistent with the 
hypothesis that the ERRs are constitutive transcriptional activa- 
tors, in vitro transcription assays demonstrated that ERRa and 
ERRy induce expression of target genes without addition of 
potential ligand to the media [11,12]. ERRa—/— mice display a 
significant decrease in body mass and ERRy— /— mice are 
perinatal lethal due to cardiac failure [13,14], phenotypes 
connected to the roles that both of these isoforms play in energy 
metabolism. 



The role of ERRs in bone and cartilage are also beginning to be 
investigated, with most data published on ERRa [15]. ERRa is 
expressed in proliferating chondrocytes in vivo and throughout 
chondrocyte differentiation in vitro [16]. In addition, it has been 
shown that ERRa is dysregulated in murine models of inflamma- 
tory arthritis [17], as well as in human osteoarthritis [18]. There is 
no data in the literature describing ERRy in cartilage, and only 
very little described on its role in bone. Results of an 
epidemiological study of ERRy polymorphisms in humans 
indicated a correlation between a subset of ERRy variants and 
elevated bone mass [19]. In vitro overexpression of ERRy causes a 
decrease in the expression of bone sialoprotein (Bsp) and 
osteocalcin (Ocn), markers of mature osteoblasts, in the MC3T3 
pre-osteoblast cell line [20]. Taken together, these results suggest 
that ERRy is a negative regulator of osteogenesis. To determine 
whether ERRy also has a biologically relevant function in 
cartilage, we have generated transgenic (Tg) mice with a collagen 
al (II) (Col2) promoter driving expression of a full length ERRy2 
(long isoform) transcript (Col2::ERRy2FL). We report here that 
overexpression of ERRy2 in a cartilage-specific manner leads to 
abnormalities in the axial and appendicular skeletons. 
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Results 

Overexpression of ERRy2 Results in Dwarfism 

To begin to investigate the putative role of ERRy in 
chondrogenesis, we first asked whether ERRy is expressed in 
cartilage. We found that ERRy is expressed in mouse cartilage at 
levels similar to ERR/} and ERfi, but approximately 50 fold less 
than ERRa, and 85 fold less than ERa (Figure 1A). To evaluate the 
consequences of cartilage-specific overexpression of ERRy, we 
generated two independent Tg mouse lines, Line 1 (#6486) and 
Line 2 (#4094), that express the longer protein isoform of ERRy2 
under the control of the collagen al (II) promoter (Co- 
12::ERRy2FL) (Figure IB). LacZ staining of 14.5dpc embryos 
revealed distinct staining in the developing craniofacial skeleton 
and both the axial and appendicular skeletons of Tg mice but not 
wild type (WT) embryos (Figure IB). Western blot analysis of 
14.5 dpc limbs demonstrated increased expression of ERRy2 
protein relative to the (3-ACTIN control and quantification 
demonstrated approximately 3-fold overexpression of ERRy2 in 
Line 1, and 5-fold overexpression in Line 2 (Figure 1C). Mice were 
born in the expected Mendelian ratios and up to at least 8 months 
of age appeared healthy, suggesting that there is no overt 
detriment due to the integration or expression of the transgene. 

Significant anomalies in all endochondral bones studied and 
measured were already apparent at birth in Tg versus WT mice 
(data not shown), but data are reported here only for postnatal day 
7 (P7). Measurement of P7 pup weight showed no significant 
difference in Line 1, but a significant reduction in Line 2 
compared to WT mice (Figure 2A). Quantification revealed a 
small but significant decrease in total body length (data not shown) 
as well as in the crown-rump length (Figure 2B) in pups of both 
lines. In P7 whole mount skeletal preparations, Tg mice also 
exhibited a shortened snout and domed head (Figure 2C) and a 
significant decrease in skull length, but not width compared to WT 
mice (Figure 2C). Analysis of specific long bones showed a 
reduction in the length of the femur and the tibia (Figure 2D) in 
both lines, as well as the humerus, ulna, radius and scapula (length 
but not width was affected) (data not shown). No significant 
differences were observed in Tg versus WT bones that form by 
intramembranous ossification, e.g., the mandible and the clavicle 
(data not shown). To rule out the possibility that the phenotypes 
we observed were due to the sex of the pups, we genotyped them 
for gender, and found that the phenotypes we describe occur in 
both male and female mice. Taken together, the data indicate that 
targeted ERRy overexpression in cartilage results in mild 
dwarfism. 

Overexpression of ERRy Impairs Chondrocyte 
Proliferation, Differentiation-maturation, Cartilage Matrix 
Production and Growth Plate Organization 

Quantification of P7 proximal humerus, distal femur, and 
proximal tibia revealed trends or significant decreases in total 
growth plate height as well as the proliferative and hypertrophic 
zones in Tg versus WT mice (Figure 3 A and B). The most 
dramatic decrease was in the proliferative zone, which displayed a 
22% height reduction in the transgenic versus WT mice, whereas 
hypertrophic zone changes were less pronounced and not 
detectable in all bones (Figure 3B). The quantification of zone 
heights in the Tg mice was complicated by 2 factors: a smaller and 
disorganized proliferative zone, and the presence of acellular 
swaths of cartilage matrix that often spanned the resting, 
proliferative and hypertrophic zones (Figure 3A). Although we 
observed 2 WT samples containing acellular swaths, such areas 
were much smaller and less pronounced than those in Tg samples. 



To determine the basis of the reduction in growth plate zone 
heights, we quantified the number of proliferating cells by 
immunostaining for the proliferation marker, Ki67 (Figure 4A). 
No significant difference was observed in resting zone chondro- 
cytes, but a 30% decrease in Ki67-positive cells was seen in the Tg 
compared to WT proliferative zone (Figure 4B). Since acellular 
masses of matrix have been reported previously in other 
genetically-modified mice, including ones with unusually wide or 
generally disorganized growth plates in which hypoxic cell death 
occurs [1], we next performed a TUNEL assay. No significant 
difference in TUNEL staining was observed in either the 
hypertrophic or proliferative zones of Tg versus WT growth 
plates (quantified for Line 2; Figure 4C). RT-qPCR also revealed 
no differences in expression of the apoptosis-associated markers, 
Bax and Bcl2, in RNA isolated from the growth plate of WT versus 
Tg mice (Figure 4D). The data suggest that the reduced length of 
the proliferative zone is a consequence of decreased chondrocyte 
proliferation, whereas that of the hypertrophic zone may be due to 
a disruption in chondrocyte differentiation, a delay in chondrocyte 
maturation, or secondary effects to the reduction in proliferation. 

To investigate the molecular basis of changes observed in the 
growth plate as a consequence of overexpression of ERRy2, we 
next used an in silico database search to screen for putative ERR 
binding sites in the region spanning 10 Kb upstream to 5 Kb 
downstream of the transcriptional start site in a variety of genes 
known to be involved in chondrogenesis. The genes that met both 
criteria are listed in Table 1 and included transcription factors 
(Sox9, Sox6, Atf4, and Rmx2), extracellular signaling molecules (Ihh, 
Bmp7, and Pthrp), hormone and growth factor receptors (Pthlr and 
Fgfr3), and cartilage matrix proteins (Col2 [Col2al], CollO 
[CollOal], and Agg). Because we observed a reduction in 
proliferation within the proliferative zone of transgenic animals, 
and because ERRy has been shown to suppress proliferation [21], 
we also assessed the expression of cell cycle regulators cyclin D 1 
[CcnDl), and cyclin dependant kinase inhibitor la and lb (Cdknla 
and Cdknlb). RT-qPCR revealed a significant increase in Pthlr 
(Figure 5A), but no differences in the transcription factors tested 
(Figure 5B). We also observed an increase in cartilage matrix 
proteins Col2 and Agg (Figure 5C) and the cell cycle regulator 
Cdknlb (p27) (Figure 5D) in Line 2 (higher overexpresser). This 
suggests that overexpression of ERRy decreases chondrocyte 
proliferation through regulation of Cdknlb, impacts growth plate 
organization, matrix synthesis as evidenced by upregulation of 
Col2 and Agg, and affects chondrocyte maturation through 
upregulation of Pthlr. 

Discussion 

We report here that overexpression of ERRy2 in chondrocytes 
results in decreased axial and appendicular skeleton size and 
disruption of growth plate height and organization. This 
phenotype was manifested already in newborn mice, indicating 
an effect on embryonic skeletal development that persisted 
postnatally, and included changes in chondrocyte proliferation, 
differentiation, and maturation, as well as matrix production, 
suggesting a role for ERRy in chondrogenesis. 

In the two independent Col2::ERRy2FL mouse lines analyzed, 
ERRy2 is overexpressed at moderate levels (3 to 5 fold higher 
levels than endogenous ERRy expression in WT littermates) and 
the resultant skeletal phenotype and concomitant changes in gene 
expression observed are relatively subtle and ERRy dose- 
dependent. It therefore seems likely that the phenotype seen in 
Col2::ERRy2FL mice is a direct consequence of modulation of 
ERRy transcriptional activity rather than from disruption of the 
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Figure 2. Skeletal analysis of P7 animals. (A) Body weight measured in P7 pups from WT (n = 1 9), Line 1 (n = 1 1 ) and Line 2 (n = 9) animals. (B) 
Crown-rump length measurement showed reduced axial skeleton. Measurements were taken from the snout to the base of the tail from WT (n = 24), 
Linel (n = 16) and Line 2 (n = 13) pups. (C) Close-up photographs of the skulls from WT and Tg animals, and resulting quantification showing reduced 
length, but not width, in Tg compared to WT animals, WT (n = 6) and Line 1 (n = 1 1 ) pups. (D) Alcian blue/Alizarin red double stain of WT and Line 2 
hindlimbs. Quantification shows reduced bone length and mineralized component of the axial skeleton in Tg compared to WT animals. Total femur 
and tibia length were measured as well as the length of the mineralized portion of the bone as demarcated by the Alizarin red staining for WT 
(n = 22), Line 1 (n = 14). Line 2 (n = 14). Graphs represent the mean ± SD *p==0.05, **p<0.01, ***p<0.001. 
doi:1 0.1 371 /journal.pone.0081 51 1 .g002 



transcriptional machinery. Nevertheless, it remains possible that 
the effects of ERRy2 overexpression are indirect, e.g., resulting 
from an imbalance in the ratio of the ERRyl/ERRy2 protein 
isoforms and sequestration of required transcription cofactors (see 
also below). Thus, further studies in ERRy knockout mice, in 
particular chondrocyte-specific ERRy knockout mice to circum- 
vent the perinatal lethality seen in global ERRy knockout mice 
([14] and Cardelli and Aubin, unpublished data), and in 
Col2::ERRyl overexpressing mouse lines, are of interest. 

The most pronounced effect observed in the Tg growth plate 
was the reduction in height of the proliferative zone, with a 22% 
reduction in proliferative zone length and a 30% decrease in the 



percentage of proliferative Tg compared to WT cells. Despite the 
fact that CmDl contains 2 putative ERREs in its regulatory region 
(Table 1 ), we did not observe any difference in CcnDl expression in 
WT versus Tg 14.5 dpc limbs. However, we observed increased 
expression of the cyclin-dependent kinase inhibitor Cdknlb in Tg 
limbs, suggesting regulation of chondrocyte proliferation through 
regulation of this kinase inhibitor, a possibility consistent with data 
showing that ERRy suppresses S-phase progression in an in vitro 
model of prostate cancer through transactivation of Cdknla and 
Cdknlb [21]. A scan of the Cdknlb regulatory region reveals 3 
putative ERRE binding motifs at positions —5077, —2625, and 
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Figure 3. Growth plate analysis of P7 pups. (A) Growth plates from WT and Line 2 femurs showing the reduced proliferative zone and the 
acellular swath observed in the transgenic animals. (B) Analysis of total, resting, proliferating, and hypertrophic growth plate heights from proximal 
humerus (PH), distal femur (DF) or proximal tibia (PT). The largest decrease in growth plate height was observed in the proliferating zones of Tg 



animals compared to WT. *p£0.05, **p; 
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;0.01, ***p<0.001 RZ, Resting Zone; PZ, Proliferative Zone; HZ, Hypertrophic Zone; AS, Acellular swath. 



— 2618, suggesting potential direct regulation, however we cannot 
rule out a novel, indirect protein interaction. 

In addition to the changes in cell cycle regulators and 
proliferation, we observed slightly reduced hypertrophic zone 
length in some bones, marked disorganization of the proliferative 
and hypertrophic zones, and acellular swaths that spanned a large 
portion of the Tg growth plate, which were not due to an increase 
in apoptosis. Moreover, Col2, Agg, and Pthlr expression were 
increased in 14.5 dpc Col2::ERRy2FL mice. Taken together, the 
data suggest that ERRy plays a role in coordinating chondrocyte 
proliferation-differentiation-matrix synthesis, but whether the 
changes are directly and causally related remain to be determined. 
Nevertheless, increased COL2, a cartilage extracellular matrix 



(ECM) protein that binds to integrins and activates signaling 
pathways essential to chondrocyte proliferation, may contribute to 
both acellular swaths of cartilage and decreased chondrocyte 
proliferation, via a negative feedback loop between abundant 
matrix production and proliferation. For example, it has been 
shown in chondrocytes from mouse ribcage that lack of P 1 integrin 
results in reduced chondrocyte motility and COL2 adhesion, as 
well as reduced CcnDl, and increased Cdknla expression, resulting 
in decreased chondrocyte proliferation [22]. Further, it has been 
shown that ERRa is involved in osteoclast migration and 
adhesion, in part through regulation of (33 integrin [23]. It is 
necessary to determine whether integrin expression and/or 
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Figure 4. Decreased proliferation, but not apoptosis, is responsible for reduction of growth plate height. (A) Hoechst labeling and 
immunofluorescence of proximal humerus growth plate from WT and Line 2 samples, showing Ki67 positive cells in the proliferating zone (PZ) (red 
boxes), and resting zone (RZ) (white box). 2° Ab is a section stained with secondary antibody only and Hoechst. Lower panel shows the proliferative 
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zone (red box) and part of the resting zone (blue box) under higher magnification, highlighting the difference in Ki67 positive cells within the 
proliferative zone of WT and Tg mice (B) Quantification of Ki67 positive cells from Line 2 proximal humerus growth plates demonstrated a clear 
reduction in the number of proliferating cells in the proliferating zone compared to the resting zone. WT (n = 6), Line 2 (n = 6); ***p£0.001 (C) 
Quantification of TUNEL positive cells shows no difference in apoptotic cells in either the proliferating or hypertrophic zones. (D) Gene expression for 
apoptotic markers, Box and Bc/2, from 14.5 dpc limbs confirms lack of increased apoptosis. 
doi:1 0.1 371 /journal.pone.0081 51 1 .g004 



chondrocyte adhesion and motility are affected in the ERRy- 
overexpressing mice. 

Whether Col2 is regulated directly by ERRy or indirecdy 
through interaction with a factor upstream of Col2 remains to be 
determined. It has been shown that overexpression of NKX3.2 
in vitro increases Col2 expression in a SOX9-independent manner, 



Table 1. Results of in silico search for putative ERRE sites in 
known regulators of growth plate chondrocytes (RZ = resting 
zone; PZ = proliferative zone; PHZ = pre-hypertrophic zone; 
HZ = hypertrophic zone). 



Position 



Gene 


Expression 


Function 


Putative ERRE (bp) 


Atf4 


PZ & PHZ 


Proliferation 


ACAAGGACA 


-7522 








TTAAGGTCA 


-6669 








ACAAGGACA 


-3214 








CAAAGGTCA 


-606 


Bmp7 


PHZ 


Proliferation & 


TAAAGGTCA 


-4858 






Differentiation 


TCTAGGTCA 


-4235 








TAAAGGTCA 


+1692 








CAAAGTTCA 


+3245 








CTAAGGTCA 


+4796 


CcnD1 


RZ, PZ, PHZ, HZ 


Proliferation 


TCGAGGTCA 


-4175 








TCTAGGTCA 


+1117 


Fgfr3 


PZ, PHZ, HZ 


Proliferation & 


GAGAGGTCA 


-6661 






Differentiation 


TGAAGGACA 


+4463 


Ihh 


PHZ 


Proliferation & 


TAAAGGTCA 


-5988 






Differentiation 


TCAAGGTCA 


-2716 








TCAAGGACA 


+1353 








CAAAGGTCA 


+4046 


Pthlr 


PZ, PHZ 


Proliferation & 


CAAAGGTCA 


-9433 






Differentiation 


CCAAGGTCA 


-2698 


Runx2 


PHZ, HZ 


Differentiation 


CAAAGTTCA 


-6895 








CAAAGGTCA 


-3156 








CAAAGTTCA 


+4145 


Sox6 


RZ, PZ, PHZ 


Proliferation & 


TCAAGGACA 


-6264 






Differentiation 


ATGAGGTCA 


-2802 








TCAAGGTCA 


-978 


Sox9 


RZ, PZ, PHZ 


Proliferation & 


CCAAGGTCA 


-8797 






Differentiation 


CCAAGGTCA 


-2966 


Pthrp 


RZ, PZ, PHZ 


Proliferation & 


TTCAGGTCA 


-2728 


Differentiation 


Co/2 


RZ, PZ, PHZ 


Cell-matrix 


TAGAGGTCA 


-8714 






communication 


CAGAGGTCA 


-7002 








GTAAGGTCA 


-2325 



Position is relative to the start site. 
doi:1 0.1 371 /journal.pone.0081 51 1 .t001 



by directly binding to a 48 bp chondrocyte-specific enhancer in 
the Col2 regulatory region [24]. Scanning the regulatory region 
using a core 'AGGTCA' sequence reveals 3 putative ERREs 
through which ERRy could directly control Col2 expression. 
Alternatively, ERRy may regulate Col2 expression through the 
recruitment of cofactors, such as PGCl-a and CBP/P300, as has 
been shown in SOX9-dependent regulation of Col2 [25,26]. 
Interestingly, it has been shown that PGCl-a is a cofactor of 
ERRy [27], and it is possible that ERRy and SOX9, together with 
PGCl-a, are part of a larger transcriptional complex to regulate 
Col2 expression. Intriguingly, one of the ERREs in the C'ol2 
regulatory region is near the SOX9 binding site, supporting the 
hypothesis of a larger transcriptional regulatory complex. 

Analyses of several genetically-engineered mouse models have 
revealed the importance of FGF [28,29] and IHH [3,30,31] 
signaling in chondrocyte proliferation and differentiation. Recent- 
ly, it has been reported that several transcription factors regulate 
the transcription of Ihh, including RUNX2 [32], ATF4 [33], and 
MSX2 [34]. The phenotype we observed in the ERRy-overex- 
pressing mice did not result in transcriptional changes in Ihh or 
Fgfr3, suggesting either a regulatory mechanism involving a 
transcriptional complex with one of the above mentioned 
transcription factors or regulation independent of IHH or FGFR3. 

It will also be important to further analyze how ERRy regulates 
chondrocyte maturation and hypertrophy, as no significant 
changes were detected in the hypertrophy markers assessed, such 
as Runx2, or its target gene CollO. While this may reflect simply the 
small content of hypertrophic cells in the 14.5 dpc samples 
utilized, it is worth noting that we also detected no differences in 
CollO expression in 17.5 dpc or P7 bone samples (data not shown). 
This suggests that ERRy may have only a small or secondary role 
in chondrocyte hypertrophy. Alternatively, the modest difference 
we observe in the Tg hypertrophic zone height may be too small to 
quantify significant differences in hypertrophic marker gene 
expression. As mentioned, the presence of the acellular region in 
the Tg animals, made measurement of the growth plate zones 
difficult, and may have influenced the observation of significant 
differences where observed. On the other hand, the increase in 
Pthlr expression seen in 14.5 dpc Tg animals may contribute to a 
reduction in the hypertrophic zone through delay in the 
differentiation process. Pthlr knockout mice exhibit shortened 
limbs characterized by decreased proliferating chondrocytes. 
However, they also display premature hypertrophy at early 
embryonic stages [35]. By contrast, transgenic mice that consti- 
tutively express Pthlr display a severe delay in the endochondral 
process, including a reduced zone of hypertrophy [36]. Although 
we did not observe any differences in hypertrophic gene expression 
in our mouse model, additional studies on the growth plate of early 
embryonic mice are needed to elucidate the basis of hypertrophic 
zone anomalies. 

In addition to the above factors, estrogen and expression of 
soluble and membrane-bound estrogen receptors (ERs and 
GPR30, respectively) have also been implicated in regulating 
growth plate chondrogenesis. Mice expressing Col2 promoter- 
driven ERot have reduced proliferation and differentiation, and 
subsequent dwarfism [6], while mice with cartilage-specific 
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Figure 5. Expression profile of putative target genes. Genes are grouped as extracellular growth factors (A), transcription factors (B), ECM 
proteins (C), and cell cycle regulators (D). Of the genes tested, we observed increased Pthh, Col2a1, Agg, and Cdknlb expression in 14.5 dpc limbs of 
Tg compared to WT animals. Expression by RT-qPCR of putative target genes was normalized to the expression of L32. WT (n = 7), Line 2 (n = 9), 
*p<0.05. 
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inactivation of ERoc exhibit prolonged longitudinal bone growth 
[5]. Taken together, the evidence suggests that not only ERa, but 
also ERRy are negative regulators of chondrocyte proliferation 
and differentiation, which appear to be opposite to the function of 
ERRot [16]. Further, ERRot is able to form heterodimers with 
ERRy [37] and ERa [38], suggesting that regulation of 
chondrocyte proliferation and differentiation may require a 
carefully controlled balance of the nuclear receptors. It remains 



to be elucidated if any interaction occurs between ERRy and 
ERRot during chondrogenesis. It is also possible that the high 
ERRy2 expression in our model may turn ERR transcriptional 
activation into repression. It was demonstrated that ERRot and 
ERRy can independently activate an ERRE-driven promoter 
reporter, but when heterodimers of ERRa/y were formed, the 
same reporter was suppressed [37]. 
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In summary, overexpression of ERRy2 in a cartilage-specific 
manner leads to dose-dependent abnormalities in the axial and 
appendicular skeletons due to alterations in Cdknlb expression and 
chondrocyte proliferation as well as differentiation-maturation- 
matrix synthesis. Work is ongoing to characterize further the 
mechanism by which ERRy exerts its actions in the developing 
growth plate. 

Materials and Methods 

Ethics Statement 

All experimental procedures were performed in accordance 
with protocols approved by the Canadian Council on Animal 
Care and the University of Toronto Faculty of Medicine and 
Pharmacy Animal Care Committee. 

Construction of pCol2a1 mERRy2 and Generation of 
Transgenic Mice 

We made a generic transgene vector containing the mouse 
Col2al promoter (gift from B. deCrombrugghe) and flanking 
intron [39] fused to a splice acceptor site [40] with stop codons in 
all reading frames, an extensive multiple cloning site (MCS), a 
Polio internal ribosome entry site (gift from P.A. Greer), a nuclear 
localized fi-galactosidase [41] and the protamine minigene pA 
[42] . During the construction of this vector, silent mutations were 
introduced to abolish restriction sites within the transgene so that 
the MCS could be expanded for future cloning of cDNAs. The 
transgene is flanked by NotI restriction sites. The longer 458 
amino acid open reading frame for ERRy2 (NM_00 1243 792) was 
PCR amplified from mouse muscle cDNA with primers containing 
restriction enzymes to facilitate cloning. All vectors were verified 
by sequencing. Transgenic lines carrying the pCol2almERRy2 
construct were generated by pronuclear injection, as described 
previously [43]. Hemizygous founders were screened for trans- 
mission of the transgene by performing PCR on DNA isolated 
from tail clips taken from the Fl generation of progeny using 
primers in the IRES (TGC TCC TTT GAA ATC TTG TGC A) 
and LacZ (AAG TTG GGT AAC GCC AGG GT) portion 
(Figure 1). To determine sex of the P0 or P7 pups, we performed 
PCR targeting the sex-determining region Y (Sry) (Table 2). 

LacZ Stain for Detection of Transgene 

Embryos of 14.5 dpc were dissected, and processed for LacZ 
detection as previously described [44]. 

Gene Expression Analysis 

Embryonic tissues and adult articular knee cartilage were 
homogenized using an Ultra Turrax T25 homogenizer, while 
postnatal tissues were manually ground with a mortar and pestle 
under liquid nitrogen. The RNA was extracted using TRIzol 
(Invitrogen), precipitated in isopropanol, and resuspended in 50- 
200 uL of DEPC dH 2 0. To remove potential contaminating 
DNA, RNA samples were subjected to DNase treatment, using an 
Ambion Turbo DNA-Free kit (Invitrogen), as per the manufac- 
turer's directions. Three ug of DNased RNA were reverse 
transcribed using Superscript II Reverse Transcriptase (Invitro- 
gen), according to the manufacturer's directions. All primers were 
designed with intron inclusion in corresponding genomic DNA, 
and are common to all potential transcript variants (Table 2). 

Western Blotting 

Dissected limbs in PBS were homogenized (Ultra Turrax T25), 
followed by lysis in RIPA buffer (50 mM Tris HC1 pH 8, 150 mM 



Table 2. Primer sequences used in gene expression analysis. 



Gene 


Upstream Sequence 


Downstream Sequence 


L32 


CACAATGTCAAGGAGCTGGAAGT 


TCTACAATGGCTTTTCGGTTCT 


Bmp7 


GCACTCAGGCAGGGAGTCGG 


ACCCAGTGGTTGCTGGTGGC 


Ihh 


TGCTGGCGCGCTTAGCAGTG 


GCAGCGGCCGAATGCTCAGA 


Pthrp 


ATTCCTACACAAGTCCCAGAG 


ACTTGCCCTTGTCATGCAGTA 


Pthh 


GACGTGGGCCAACTACAGCG 


GTGCAGTGCAGCCGCCTAAA 


Fgfr3 


TAGCGGCCGCCAGTCTCCAC 


ACGCAGGCCGGGACTACCAT 


Sox9 


AATGCTATCTTCAAGGCGCTG 


GGACCCTGAGATTGCCCAG 


Sox6 


ACAACCACAGACAGATTGAGCAGC 


TGCCCCTGCCGAGTTTGGTG 


Atf4 


ATGGCGTATTAGAGGCAGCA 


GATTTCGTGAAGAGCGCCAT 


Runx2 


TGTTCTCTGATCGCCTCAGTG 


CCTGGGATCTGTAATCTGACTCT 


Col2 


ACTGGTGGAGCAGCAAGAGC 


TCTGGACGTTAGCGGTGTTG 


Col 10 


AACGGTACCAAACGCCCAC 


CTTTGTTCTCCTCTTACTGGAATCCC 


Agg 


GCGTGAGCATCCCTCAACCATC 


GGCAGTGGTCACAGGATGCATG 


CcDI 


CCTGTGCGCCCTCCGTATCT 


TCATGGCCAGCGGGAAGACC 


Cdknla 


CAGACCAGCCTGACAGATTTCTA 


GAGGGCTAAGGCCGAAGATG 


Cdknlb 


GTTTCAGACGGTTCCCCGAA 


TCTTAATTCGGAGCTGTTTACGTC 


ERRu. 


TCGAGAGATAGTGGTCACCATCAG 


CTTCCATCCACACACTCTGCAG 


ERRfi 


TGAGATCACCAAACGGAGGC 


GAACTCGGTCAAGGCGCA 


ERRy 


TGTGACTTGGCTGACCGAGA 


TGGAGGAGGCTCATCTGGTCT 


ERa 


GGCTGCGCAAGTGTTACGAA 


CATTTCGGCCTTCCAAGTCAT 


ERji 


TTGGTGTGAAGCAAGATCACTAGAA 


GACTAGTAACAGGGCTGGCACAA 


Sry 


GAGAGCATGGAGGGCCAT 


CCACTCCTCTGTGACACT 



doi:1 0.1 371/joumal.pone.0081 51 1 .t002 



NaCl, 1% NP-40, 0.5% sodium deoxycholate and 0.1% sodium 
dodecyl sulphate) with added protease inhibitors. Protein samples 
were quantified using the Bio-Rad DC Protein Assay kit, following 
the manufacturer's instructions. Thirty u.g of each sample was run 
in a 10% SDS-PAGE gel, transferred to polyvinylidene fluoride 
(PVDF) membrane, followed by blocking in 5% milk-TBS-T for 
30 minutes at room temperature. Immunodetection was carried 
out using a rabbit polyclonal anti-ERRy antibody (H38x, Santa 
Cruz Biotechnology Inc.) diluted 1:5000 in blocking buffer, or 
rabbit anti-P-ACTIN antibody diluted to 1 :2000 (Sigma). This was 
followed by a one hour incubation with a goat anti-rabbit IgG, 
conjugated to horse radish peroxidase (HRP; Santa Cruz 
Biotechnology), diluted 1:5000-1:8000 in blocking buffer. The 
HRP was visualized using the Amersham ECL Western Blotting 
Detection kit (GE Healthcare), as per the manufacturer's 
instructions. The autoradiographic films from both ERRy and 
(3-ACTIN detection were scanned and the density and size of the 
bands were quantified using Image Lab software (Bio-Rad). The 
ERRy band was normalized to the fi-ACTIN band to assess 
proportionate protein levels. The 29 amino acid difference 
between the ERRyl and ERRy2 protein isoforms is not clearly 
obvious in the gels that we run and we attribute the increased 
expression we see to the ERRy2 isoform. 

Whole Mount Skeletal Staining 

P0 and P7 animals were dissected, eviscerated, and fixed in 95% 
ethanol overnight or up to two weeks, and then processed for 
whole mount skeletal staining as previously described [45]. When 
samples were fully cleared, skeletons were dissected and photo- 
graphed in a Petri dish containing 100% glycerol, using a Nikon 
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Coolpix P5100 digital camera affixed to a dissecting scope. The 
images were then quantified in Image J by taking linear 
measurements of individual skeletal components. 

Histological and Immunofluorescence Analysis 

The left limbs taken from each skeleton analyzed by whole 
mount staining were fixed in 4% paraformaldehyde for 24 hours, 
transferred to PBS for 2-3 days and then decalcified (10% EDTA, 
0.1 M Tris pH 7.4) for 15 days. Once decalcified, limbs were 
serially dehydrated with one 20 minute wash in 30% ethanol, 
three 20 minute washes in 50% ethanol and a final wash and 
storage in 70% ethanol. Samples were paraffin-embedded, 
sectioned (5 um), and stained (hematoxylin and eosin, and 
Safranin-O) at the Toronto Centre for Phenogenomics (TCP). 
The Safranin-O-stained slides were photographed under a Nikon 
Eclipse TS100 fluorescence microscope, at magnifications of 40X 
and 100X. In the 40X micrographs, measurements were taken 
from the absolute edge of Safranin-O staining before the articular 
surface, down to the last hypertrophic chondrocyte before the 
primary spongiosa. Three such measurements were taken at 
separate axes along the bone, with the average representing the 
total growth plate height. Using the 100X micrographs, measure- 
ments were taken from the first flattened chondrocyte, which 
appeared in a distinct column of three or more cells, down to the 
last hypertrophic chondrocyte, and then the subset of that distance 
that contained only hypertrophic chondrocytes was taken. This 
was repeated at five points along the width of each growth plate, 
and the averages were used to calculate the specific heights of the 
resting, proliferative and hypertrophic zones. For each animal, 
growth plates were assessed from the proximal humerus, distal 
femur and proximal tibia. 

To immunodetect Ki67, a common proliferation marker [42], 
sections were deparaffinized and rehydrated in ethanol washes, 
followed by antigen retrieval in boiling citrate buffer ( 1 0 mil citric 
acid, 0.05% Tween-20, pH 6) for up to 20 minutes. The slides 
were blocked in normal goat serum (Invitrogen) for 30 minutes at 
room temperature, washed, incubated with rabbit polyclonal anti- 
Ki67 antibody (diluted 1:25 in blocking buffer) for 1 hour at room 
temperature, washed, then incubated with goat anti-rabbit Alexa- 
594-conjugated secondary antibody (Invitrogen, diluted 1:50 in 
blocking buffer) for 30 minutes at room temperature. The samples 
were counterstained with Hoechst nuclear stain. 

Growth plate sections of the proximal humerus were viewed 
and imaged using a Bioquant Osteo imager with Photofluor II 
fluorescence excitation light, a triple Chroma filter, and Bioquant 
Osteo 2012. Ki67 positive and negative cells within the 
proliferative and resting zones were counted in Image J. By 
calculating the percentage of total cells (defined by the Hoechst 
stain) that stained positively for Ki67, the mitotic activity in the 
proliferative zone and resting zone was measured. 

TUNEL Assay 

Growth plate sections of the proximal humerus were processed 
as described above, before using the FragEL DNA fragmentation 
detection kit (Calbiochem), as per the manufacturers instructions, 
and counterstained with methyl green. Sections were imaged at 
16X magnification, using a Nikon Eclipse TS100 fluorescence 
microscope. Images were analyzed using ImageJ. Cells within the 
proliferative and hypertrophic zones were counted for TUNEL 
analysis and normalized to methyl green stained nuclei. 



Expression Analysis of Putative Target Genes 

A list of putative target genes was constructed by in silico analysis 
of the regulatory regions of factors known to be involved in the 
chondrocyte differentiation process and that contained at least one 
ERR binding site (ERRE) in the region from - 10 Kb to +5 Kb of 
their annotated transcription start site. The search used the ERRE 
consensus sequence TCAAGGTCA and 20 additional sequence 
variants that had been identified in the literature [46] . The search 
was performed using the freely available Transcriptional Regula- 
tory Element Database (TRED) (http://rulai.cshl.edu/cgi-bin/ 
TRED/ tred.cgi?process = home), which retrieves genomic se- 
quences from the current Ensembl build of the mouse genome. 
The results from this search are shown in Table 1 , along with the 
primers used to look for differences in gene expression in Table 2. 
The primers were chosen to pick up all of the known transcript 
variants and include at least one intron in the corresponding 
genomic DNA. 

The reactions were performed in triplicate on a 96-well plate in 
a BioRad MyIQ_ iCycler, for 50 cycles with an annealing 
temperature of 59°C. The amplification data was uploaded into 
the PCR miner program (http://www.ewindup.info/miner/ 
version2/) to obtain the Ct and reaction efficiency values. The 
relative expression levels of the target gene were normalized to the 
L32 internal control expression. 

Statistical Analysis 

All data were analyzed using Graphpad Prism 4.0 software. 
When three data sets were analysed, ANOVA was used first to 
determine significance, followed by Student's t-Test. All the graphs 
are plotted as the mean ± standard deviation and the p values 
listed are for the comparison to the WT values. Graphs were 
constructed using Microsoft Excel 2003 software. 

Mouse Gene Nomenclature 

We followed the mouse nomenclature guide as stated on the 
Mouse Genome Informatics web page (http:/ /www.informatics. 
jax.org/mgihome/nomen/short_gene.shtml). Thus, mouse genes 
are written with first letter capitalized followed by small letters, all 
italicized, while proteins are written all capitalized and without 
italics. 
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